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ABSTRACT 


The  amplification  of  P wave  amplitudes  by  the  receiver  crustal  structure 
was  modelled  for  34  LRSM  stations  by  using  crustal  structures  derived  from 
geological  and  geophysical  information  in  the  literature.  Crustal  amplifi- 
cation was  found  to  be  linearly  related  to  the  acoustical  impedance  of  near 
surface  materials  at  each  site.  Correlation  of  the  crustal  amplification, 
expressed  in  magnitude  units,  with  the  magnitude  residuals  of  Booth,  Marshall 
and  Young  (1974)  is  statistically  not  significant  for  the  complete  set  of 
stations  investigated  by  the  above  authors.  This  indicates  that  crustal 
amplification  alone  cannot  explain  the  observed  magnitude  residuals.  If  two 
subsets  of  stations  in  western  and  eastern  United  States,  respectively,  are 
taken  a linear  trend  between  the  computed  crustal  amplification  and  the 
magnitude  residuals  can  be  seen  within  each  set,  but  there  is  a separation 
of  about  .3  magnitude  units  between  the  two  groups.  The  most  likely  explana- 
tion of  the  separation  is  anelastic  attenuation  in  the  mantle  under  the 
western  United  States.  About  75%  of  variance  in  the  magnitude  residrals,  for 
the  US  stations  used,  can  be  correlated  with  anelasticity  and  crustal 
amplification  as  known  at  present.  The  residual  variance  must  be  due  to 
uncertainties  in  the  above  factors  and  other  causes. 
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Stiuilcs  of  teloso Ism Ic  nuignitiiclo  mo.'isurt'nients  have  sliown  that  systematic 
il  1 f fereni'es  exist  between  magnitudes  observed  at  various  stations  wiiich  can- 
not be  associated  with  regional  source  effects  or  radiation  patterns.  Tlie 
anomalies  must,  therefore,  he  associated  with  the  etfects  of  the  crust  and/or 
the  upper  mantle  under  each  station.  Booth,  Marsltall  and  Young  (1974)  have 
published  a set  of  station  magnitude  corrections,  mostly  for  LRSM  stations  in 
North  America.  The  purpose  of  this  study  is  to  evaluate  the  contribution  of 
crustal  structure  to  the  station  anom.-illes.  The  structure  of  the  North  Ameri- 
can continent  lias  been  studied  quite  extensively',  the  crustal  structure  is 
reasonably  well  known  at  many  hRSM  stations.  At  a few  stations,  however,  the 
striK-tiire  is  not  known  and  at  most  stations  the  fine  structure  is  not  well 
known.  In  such  c.-ises  a rough  model  can  be  obtained  by  using  the  regional 
crustal  thicknesses  and  surface  geological  Information. 

We  have  compiled  a set  of  estimated  crustal  models  for  the  majority  of 
the  stations  analyztxl  by  Booth,  et  al,  (1974).  In  this  report  we  shall  com- 
pare tlie  crustal  magnification  computed  for  each  station,  based  on  tlie  assumed 
models,  witli  tlie  observed  magnitude  anomalies.  'I’he  location  of  stations  within 
tlie  continental  II. .S.  and  the  short-period  residuals  are  shown  In  Figure  1. 
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Figure  1.  Short-period  magnitude  residuals  in  the  United  States  determined 
by  Booth,  Marshall  and  Young,  1974. 
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CRUSTAL  MODELS 


Since  we  only  want  to  consider  the  crust  we  terminate  our  modc'ls  with  a 
lialfspace  below  the  Moho,  The  properties  of  the  local  material  below  the 
Moho  are  determined  from  refraction  studies  at  stations  where  major  refrac- 
tion svirveys  are  available.  At  some  stations  the  (;eotech  LRSM  site  reports 
were  used  to  olualn  near-surface  Information.  These  reports  are  based  on 
flrstliand  observations  by  geologists,  together,  usually,  wltl\  perusal  of 
geological  maps  of  the  area.  Especially  Important  parameters  for  computing 
the  crustal  amplification  are  the  velocities  and  thicknesses  of  the  uncon- 
solidated sediments  near  the  surface.  The  lithology  by  Itself  Is  not  directly 
appreciable  in  this  study.  Althougli  many  previous  studies  stressed  the  im- 
portance of  tile  nature  of  the  rocks  at  the  surface  e.g.,  Evernden  and  Clark 
(l^>70),  ucli  presentations  can  be  very  misleading.  The  velocity  of  sedimen- 
tary rocks  is  more  important  than  their  composition.  Some  older  sedimentary 
rocks  are  extremely  compact  and  l\ave  high  seismic  velocities,  wliict\  are  com- 
parable to  those  of  granite  and  produce  amplltvides  comparable  to  granite.  Cn 
the  other  hand,  thick  low  velocity  sediments  can  considerably  Increase  the 
amplitude  of  seismic  waves  at  the  surface. 

Crustal  models  for  the  set  of  LRSM  stations  used  are  given  In  Table  1. 

As  stated  above,  these  are  of  varying  quality  and  may  he  improved  as  more 
data  becomes  available.  Besides  these,  the  Ceoteeh  LRSM  site  reports  were 
always  considered  In  constructing  the  near  surface  part  of  eacl)  model.  in- 
formation in  these  reports  varies  from  quite  detailed  to  sketchy.  A consid- 
erable amount  of  Information  about  the  noar-surface  structure  has  been  gen- 
erated by  exploration  for  oil,  hut  such  Information  Is  not  easily  available. 
There  Is  undoubtedly  more  Information  In  the  general  geological  literature, 
hut  such  information  Is  very  hard  to  ctillect,  and  the  geological  literature 
usually  does  net  give  tljc  relevant  parameters,  the  seismic  velocities  and 
densities,  necessary  to  construct  an  accurate  geopltynlcal  model.  Tlte  set  of 
mi>dels  given  in  this  report  represents  what  can  he  collecle»l  with  a modest 
effort  from  tlu'  literature.  Quality  estimates  have  been  made  for  each  sec- 

Evernden  and  Clark,  l‘)70.  .Study  of  teleselsmlc  P.  11.  amplitude  data,  IMiys. 

Earth  Planet  InterJU>r_s,  24-31. 
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tion.  If  hard,  high-velocity  rock  is  at  the  surface,  the  section  should  be 
excellent.  Soft  rock  sites  can  be  excellent  if  detailed  well-logs  are 
available.  If  larger  thicknesses  of  the  surface  soft  layer  must  be  guessed, 
then  the  section  quality  is  poor. 


TABLE  1 


Crustal  Models 

for  the  LRSM 

Stations 

in  t ))  1 s 

Study 

Section  Quality  E 

- Excellent, 

GF  - Good 

to  Fair 

, P-  Poor 

STATION 

L*y«r  Thtcknesii 

Compreititlonfll 

Ve  tuc  icy 

Shear 

Ve lor ity 

SOURCE 

(QUALITY) 

d 

a 

8 ■ 

p 

REFERENCES 

AD-IS 

6.00 

5.50 

2.80 

2.70 

Jacobs  and 

or 

21.00 

6.40 

3.60 

2,90 

Hamada  (1972) 

oo 

8.10 

4.60 

3.30 

AR-WS 

6.00 

5.64 

3.47 

2.70 

Bruno  and 

E 

10.50 

6.15 

3.64 

2.80 

Dorman  (1963) 

18.70 

6.60 

3.85 

2.85 

OO 

8.10 

4.72 

3.30 

AX-AL 

13.00 

6.00 

3.46 

2.70 

Antoine  and 

GF 

22.00 

6.70 

3,89 

2,90 

Ewing  (1963) 

OO 

8.10 

4.67 

3.30 

Woolard  (1959) 

BE-FL 

1.00 

2.60 

1.40 

1.80 

Antoine  and 

GF 

0.50 

3.60 

2.00 

2.00 

Ewing  (1963) 

15.00 

6.20 

3.60 

2.80 

Woebor  and 

17.00 

6.70 

3,87 

2.80 

Pennhollow  (1975) 

00 

8,10 

4.60 

3.30 

BL-WV 

1.89 

6.00 

3.46 

2,70 

Woolard  (1959) 

GF 

25.10 

6.15 

3.55 

2.82 

9.00 

6.58 

3.80 

2.92 

uo 

8.14 

4.76 

3.53 

BR-PA 

1.89 

6.00 

3.46 

2.70 

Woolard  (1959) 

GF 

25.10 

6.15 

3.55 

2.82 

9.00 

6.58 

3.80 

2.92 

OO 

8.14 

4.76 

3.53 

CP-CL 

4.90 

5.80 

3.38 

2.60 

llealy  (1973) 

E 

11.10 

6.00 

3.50 

2.70 

Pakiser  and 

11.10 

7,10 

4.00 

2.80 

Steinhart  (1964) 

OO 

7.90 

4.50 

3.30 

DH-NY 

1.89 

6.00 

3.46 

2.70 

Woolard  (1959) 

GF 

25.10 

6.15 

3.55 

2.82 

9.00 

6.58 

3.80 

2.92 

< 

OO 

8.14 

4.76 

3.53 

UR-CO 

2.00 

6.00 

3.42 

2.70 

Keller,  Smith 

E 

7.00 

6.10 

3.52 

2.70 

and  Braile  (1975) 

15.00 

6.40 

3.69 

2,80 

Keller,  Smith, 

16.00 

6.80 

3.92 

2.90 

Braile,  Heaney, 

OO 

7.80 

4.50 

3.30 

Shurbet  (1976) 
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TABLE  T (Continued) 


Crustal  Models 

for  the  LRSM 

Stations 

In 

this 

Study 

Section  Quality  E - 

Excellent , 

GF  - Good 

to 

Fair, 

P - Poor 

STATION 

(QUALITY) 

L«y«r  Thickn««* 

Coaprcoatonal 
V«luc Ity 

Shi'*r 

VclQcUy 

SOURCE 

REFERENCES 

d 

a 

8 

p 

KK-CO 

0.60 

2.90 

1.55 

1.80 

Jackson,  Steward 

GF 

0.40 

4.80 

2.64 

2.40 

and  Pakisor  (1963) 

1.30 

5.20 

3.00 

2.60 

9.90 

5.80 

3.34 

2.70 

15.50 

6.10 

3.52 

2.80 

20.00 

6.70 

3.89 

2.90 

07 

8.00 

4.62 

3.30 

FM-UT 

1.70 

3.40 

2.00 

2.30 

Keller,  Siiiitli,4 

GF 

6.70 

6.00 

3.50 

2.80 

Braile  (1975) 

6.30 

5.50 

2.90 

2.70 

Keller,  Smith, 

10.00 

6.50 

3.50 

2.80 

Braile,  Heaney, 

CP 

7.40 

4.00 

3.10 

& Sliurbet  ( 1976) 

GV-TX 

1.00 

3.05 

1.63 

1.80 

Cram  (1961) 

E 

1.00 

4.11 

2.26 

2.20 

TryuRvason  S. 

1.00 

4.42 

2.55 

2.30 

Qualls  (1967) 

6.30 

5.63 

3.24 

2.70 

15.90 

6.10 

3.52 

2.80 

16.80 

6.92 

4.00 

2.90 

tio 

8.20 

4.73 

3.40 

HL-ID 

8.00 

5.20 

2.96 

2.50 

Hill  4 Paklser 

GF 

37.00 

6.70 

3.82 

2.80 

(1966) 

CP 

7.90 

4.50 

3.20 

HN-ME 

10.00 

5.90 

3.36 

2.70 

Dainty,  Keen, 

GF 

20.00 

6.35 

3.62 

2.72 

& Blanchard  (1966) 

25.00 

7.35 

4.19 

2.83 

00 

8.10 

4.62 

3.30 

HW-IS 

3.50 

3.70 

2.11 

1.94 

Hill  (1969) 

E 

6.50 

6.00 

3.42 

2.70 

5.00 

7.20 

4.67 

3.40 

oo 

8.20 

4.67 

3.40 

JE-LA 

1.50 

2.70 

1.44 

2.00 

Antoine  and 

E 

6.60 

3.60 

2.00 

2.30 

Ewing  (1963) 

13.00 

5.80 

3.37 

2.70 

Woeber  and 

15.00 

6.70 

3.90 

2.80 

Pennhollow  (1975) 

1 

OO 

8.10 

4.80 

3.30 

1 

JR-AZ 

0.36 

2.50 

1.33 

1.80 

vU'hnson  (1965) 

GF 

19.00 

6.00 

3.50 

2.70 

12.00 

6.70 

3.90 

2.80 

CD 

7.90 

4.50 

3.30 
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TABLE  T (Continued) 


Crustal  Models 

for  the  LRSM 

Stations 

in  this 

Section  Quality  E 

- Excellent, 

CF  - Cood 

to  Fair 

STATION 

L«y«r  Thlckti«»« 

CuMpresslunal 

V«loc 1 ty 

Shi'«r 

Vrlorlly 

ilrniilty 

( QUALITY) 

d 

a 

8 

p 

KC-MO 

1.9A 

5.00 

2.88 

2.60 

CF 

9.88 

6.06 

3.50 

2.70 

12.35 

6.27 

3.62 

2.80 

18. 3A 

6.71 

3.87 

2.90 

W 

8.09 

4.70 

3.30 

KN-UT 

2.00 

3.15 

1.81 

2.00 

CF 

7.00 

6. 10 

3.52 

2.60 

15.00 

6.40 

3.69 

2.70 

16.00 

6.80 

3.92 

2.80 

at 

7.80 

4.50 

3.20 

LC-NM 

1.20 

3.00 

1.60 

2.20 

P 

5.00 

5.80 

1.20 

2.70 

21.00 

6.50 

3.80 

2.90 

ot* 

7.90 

4.45 

3.30 

LV-Ij\ 

1.20 

2.70 

1.44 

2.00 

CF 

6.00 

3.60 

2.00 

2.30 

13.90 

5.80 

1.37 

2.70 

15.00 

6.70 

3.90 

2.80 

at 

8.10 

4.80 

3 . 30 

MN-NV 

2.00 

5.80 

3.34 

2.50 

CF 

27.00 

6.15 

3.54 

2.70 

10.00 

7.10 

4.09 

2.90 

IV 

7.80 

4.50 

3.30 

MO- ID 

0.55 

1.50 

0.75 

1.80 

CF 

0.43 

4.00 

2.20 

2.50 

7.50 

5.20 

2.97 

2.80 

37.00 

6.70 

3.82 

2.80 

7.90 

4.50 

3.20 

MV-CL 

1.00 

3.00 

1.60 

2.20 

P 

15.00 

6.00 

3.50 

2.70 

20.00 

6.70 

3.90 

2.90 

CO 

7.80 

4.40 

3.30 

NC-WS 

6.00 

5.64 

3.47 

2.70 

CF 

10.50 

6.15 

3.64 

2.80 

18.70 

6.60 

3.85 

2.85 

00 

8.10 

4.72 

3.30 

Study 
P - Poor 

SOURCE 

REFERENCES 

Stewart  (1968) 


Keller,  Smith 
& Braile  (1975) 
Keller,  Smith, 
Braile,  Heaney, 
Shurbet  (1976) 

Toppozada  and 
Sanford  (1976) 


Antoine  and 
Ewing  (1963) 


Hill  & Pakiser 
(1966) 


Hill  & Pakiser 
(1966) 

(installed  in 
a missile  silo 
on  silt.) 

Pakiser  & Zietz 
(1965) 


brune  & Dorman 
(1963) 
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TART,F.  T (Contlmied) 


Crustal  Models 

for  the  I.RSM 

Stations 

In  this 

Study 

Section  Qiialltv  F 

- Excellent, 

GF  - Good 

to  Fair 

, F - Poor 

STATION 

(QUALITY) 

Thickness 

d 

CtfMpi ion«l 

V* loi  tty 

L) 

VvliK  ttV 

8 

unity 

SOURCE 

REFERENCES 

NF-NT 

1 . 00 

4.  13 

2 . SO 

2.40 

Brune  i»  l)orm^m 

c:k 

9 . 00 

h . 00 

l.SO 

2 . 70 

(1983) 

2h.00 

7.  30 

4.  IS 

2.90 

8.20 

4,70 

3 . 30 

FG-BC 

2 . 00 

5 . 00 

2.90 

2.80 

White  iS  Sava>;e 

F 

14.00 

8.10 

3.50 

2.70 

(1985) 

14.00 

8.  70 

3.84 

2.80 

w 

8.00 

4.80 

3.30 

RK-ON 

6 . 00 

S.84 

3.47 

2.70 

Bruiie  \ Dorman 

K 

10.  SO 

8.  IS 

3 . 84 

2.80 

(1981) 

18.70 

8.80 

1.8S 

2.8S 

vV 

8.  10 

4.72 

3.  10 

SJ-TX 

1 . 00 

2.13 

1.13 

1.80 

Woeher 

E 

S . U) 

3.  IS 

1.82 

2 . 00 

Fenuhollow  (197S) 

1.70 

3 . 80 

2.07 

2.40 

11.00 

8 . 00 

1.48 

2.70 

1 1 . 00 

8.70 

3.87 

2 . 90 

8 . 00 

4.82 

3 . 30 

SV-QB 

E 

1.13 

S.80 

3 . 38 

2.80 

Brune  & Donn.in 

20.00 

8.20 

3 . 80 

2.70 

(19b3) 

18.00 

8.70 

3 . 90 

2.90 

8 . 00 

4.80 

1 . 30 

SM-MA 

4.70 

S.OO 

2. 85 

2 . SO 

Asada  N Aldrich 

GF 

12,00 

8.20 

3.51 

2.70 

(19h8) 

39.00 

8.70 

3.82 

2.80 

8 . 00 

4.S8 

3 . 30 

TK-Gl, 

2.00 

3 . 00 

1.81 

2 . 00 

Mealy  (1971) 

F 

14.72 

8.10 

3.52 

2.70 

Faklser  & Zletz 

9.40 

7 . 00 

4.04 

2.90 

(1985) 

8.  10 

4.87 

3 . 30 

WH-YK 

IS.  00 

8.20 

1.54 

2.75 

LeBlanc  i 

F 

rt* 

8 20 

4. 75 

3 . 30 

Wetmiller  (1974) 

Wl-NV 

1 . 00 

4.  SO 

2,57 

2 . 10 

GK 

18.00 

8.00 

3.42 

2.70 

12.00 

8.70 

3.82 

2.80 

rt' 

7 . 90 

4.  SO 

3 . 20 
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CXWIITATIONS 


Till.'  I'nistiil  amplification  Is  a fvinction  of  not  only  the  tmitorlal  property 
distribution  but  also  ol  the  angle  of  Incidence.  We  evaluated  the  ground 
amplification  for  t)iree  different  pulse  shapes  using  an  angle  of  Incidence 
of  25“  below  the  Moho,  wtilcli  Is  fairly  typical  for  teleselsmlc  arrivals.  Tlu* 
first  pulse  Is  citosen  to  represent  the  wide-hand  signal  generated  by  a 5 kt 
nuclear  explosion  In  granite  as  seen  tltrougli  the  standard  LRSM  short-period 
Instrument  (SP5),  the  second  Is  a 50  kt  explosion  seen  througl»  the  same  in- 
strument (Sl’5())  and  the  third  Is  the  50  kt  explosion  seen  through  the  l.RSM 
long-period  Instrument  (1.1’50). 

At  first  we  performed  exploratory  calculations  and  found  that  the  ampli- 
fication is  relatively  Insensitive  to  the  moderate  variations  of  the  angles 
of  Incidence  expected  for  the  teleselsmlc  distance  range.  As  we  shall  see  In 
detail  In  tills  report  the  amplification  was  also  found  to  be  almost  as  well 
predicted  by  tlie  elastic  properties  of  materials  near  the  surface  as  by  tlie 
detailed  variations  of  tlie  crustal  frequenev  response  with  frequency.  The 
dominant  period  of  tlie  transmitted  pulses  was  also  found,  by  visual  Inspec- 
tion, to  be  Insensitive  to  crustal  structure.  Anticipating  results  to  come, 
comparison  of  the  amplitudes  of  station  pairs  wliose  structures  diffi'reil  by 
5-lOZ  e.g.,  ADIS-AIOM.,  ARWS-AXAl,,  l.ONM-MVCL,  OUNY-DRCO,  NPNT-WINV;  reveals 
iTkignitude  differences  ranging  from  .01  to  .05.  Thus,  we  feel  that  these 
crustal  models  are  about  as  accurate  as  needed  for  tlie  purposes  of  estimating 
magnit  lilies  since  most  of  the  structures  should  be  good  to  lOX.  Variations 
due  to  focussing  should  Introduce  more  variance  than  sliould  Inaccruate  plane- 
layer  parameters. 

The  algorltlim  to  compute  crustal  amplification  is  the  well-known  Haskell 
matrix  metliod  (Haskell,  ldh2)  for  laterally  homogeneous  elastic  layered  media. 
No  anelastlc  effects  In  tlie  crust  were  considered.  We  assume  a dilatational 
strain  amplitude  (W  » — + wliere  u,  v,  w are  the  components  of 

displacement)  of  unity.  In  Haskell's  notation  A"  =•  I,  for  the  incident  wave 
at  the  base  of  our  model.  This  amounts  only  to  a convenient  normalization  ol 
the  Incident  amplitude.  To  equalize  the  elastic  energy  1 lux  at  the  base  ol 

i 1/2 

the  moilel  the  results  of  computation  have  to  be  divided  bv  (p  a ) where 

' n n 

Haskell,  1H62.  Crustal  reflection  of  P and  SV  waves,  C.eopliys . Res. , b7 , 
4751-47f)7. 
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and  are  the  density  and  P wave  celocity  in  the  sub-Moho  halfspace  ter- 
minating the  model.  This  correction  factor  can  easily  be  derived  from 
Haskell's  equations  as  given  by  Bath  (1968,  p.  282-283)  and  Ewing,  Jardetzky, 
and  Press  (1957,  p.  28).  Figure  la  shows  the  location  of  LRSM  stations  used 
within  North  America. 

Figures  2 to  30  show  the  crustal  frequency  responses,  plotted  as  power 
spectra,  of  the  impulse  response  of  the  crust;  time  domain  crustal  impulse  re- 
sponses, and  final  pulse  shapes  for  each  model  and  assumed  pulse  type.  The 
LRSM  station  in  question  is  Indicated  at  the  top  of  each  figure.  Below  this, 
the  crustal  response  is  shown  as  a function  of  frequency  between  0-6  Hz.  The 
short-period  and  long-period  impulse  responses,  and  synthetic  seismograms  are 
shown  below  these  with  the  appropriate  time  scales.  All  plots  are  normalized 
to  the  same  amplitude  and  the  normalizing  factors  are  written  beside  each  plot. 
Signals  which  seem  to  be  precursors  should,  in  reality,  follow  the  first  main 
pulse  and  were  created  by  the  "wraparound"  affect  of  the  Finite  Fourier  trans- 
form. The  effect  of  these  on  the  logarithms  of  the  main  pulse  is  negligible  i 

since  they  are,  in  all  cases,  very  small.  The  maximum  amplitude  is  picked  by 
the  computer  and  printed  on  the  plot.  Only  the  relative  changes  of  amplitude  ^ 

for  a fixed  pulse  type  are  significant,  not  the  absolute  value. 

The  receiver  crust  frequency  responses  show  rapid  fluctuation  with  fre- 
quency, but  averages  over  wider  frequency  ranges  are  approximately  constant 
in  spite  of  the  fact  that  the  power  is  plotted  which  tends  to  emphasize  the 
fluctuations  on  the  linear  scales  employed.  In  the  overwhelming  majority  of 
the  cases  presented,  the  effect  of  crustal  responses  in  the  estimation  of  Q 
by  fitting  straight  lines  to  spectral  ratios  would  seem  to  be  negligible, 
since  there  is  no  clear  trend  with  increasing  frequency.  The  few  exceptions 
(such  as  FKCO)  serve  to  demonstrate  that  crustal  responses  can,  in  some  cases, 
contribute  to  the  scatter  in  Q measurements. 

— — — t 

Bath,  Marcus  (1968).  Mathematical  aspects  of  seismology,  Elsevier  Publishing 

Company,  Amsterdam. 

Ewing,  M.  W.,  Jardetzky,  W.  S.,  and  F.  Press  (1957).  Elastic  waves  in  layer- 
ed media.  McGraw-Hill  Book  Co.,  New  York,  NY 
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Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  Impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  A1)-IS.  The 
numbers  to  the  right  of  abscissa  Indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  3.  Crustal  frequency  response  function,  shov«ni  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  AX-AL.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Fij;urt'  4.  Crustal  frequom-y  response  function,  sl>owi\  as  linear  pi'wer  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  l,RSM  station  HK-Kl,,  The 
numbers  to  the  rlsht  of  absciss.!  Indicate  tlie  maximum  power  i>r 
amplitude  as  appropriate. 
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Figure  5.  Crustal  frequency  response  function,  shovm  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  CP-CL.  The 
numbers  to  the  right  of  abscissa  Indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Crustal  frcciuoncy  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  tlie  l.RSM  station  B1*-WV,  UR-l’A 
and  Dll-NY.  Tlie  numbers  to  the  right  of  abscissa  indicate  the  max- 
imum pi>wer  or  amplitude  as  appropriate. 
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CRUSTAL  RESPONSE 
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Klsiiro  1,  CriiNtal  fri'<]iu'ncy  rospDiiso  liini'tlon,  slu'wn  jin  linoiir  pnwi'r  spoi-t  mm 
ot  tho  roci'lvor  cmstJil  Impiilso  rosponso,  rocoi vor  i-mst  Impuls*' 
ri'sponsos,  synthetic  si'lsmourjims  lor  tho  I.KSM  stalioii  l>K-(’0.  Tho 

niinibors  to  tlio  right  of  iibsolsNa  liiiliiMto  tho  maximum  pow«>i  oi 
jimplltiiilo  as  jippropr late. 
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Klgiiro  8.  OniHtal  friMiuonry  response  tnnetlon,  shown  as  linear  pi>wer  speetrnm 
of  tl>e  receiver  crustal  impulse  response,  receiver  crust  impulst' 
responses,  synthetic  seismoprams  tor  the  LKSM  statli>n  I'K-tXi.  Tlu' 
numhers  to  the  rlpht  ot  abscissa  Indicate  tlu'  maximum  power  or 
amplitude  as  appropriate. 
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Figure  9.  Crustal  frequency  response  funetlon,  shown  as  linear  power  speetrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  Impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  KM-H'l'.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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10,  C.nistal  friMiuoncy  rosponso  f'uni'tion,  slu'wn  as  linear  pnwor  spi'i't  rum 
of  the  rei'i'ivor  crustal  impulse  response,  reci'iver  crust  impulse 
responses,  synthetic  seisnu'nrams  for  the  l.KSM  station  OV-i'X.  Tlie 
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Klgure  11.  Crustal  froquency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  Impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  HL-ll).  The 
numbers  to  the  right  of  abscissa  Indicate  tl»e  maximum  power  or 
amplitude  as  appropriate. 
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Figure  12.  Crustal  frequency  response  function,  sliown  as  linear  power  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  Impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  HN-MK.  The 
numbers  to  tlie  right  of  abscissa  lndic;ite  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  13. 


Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  HW-IS.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  14.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  AR-WS,  NC-WS, 
PM-WY  and  RK-ON.  The  numbers  to  the  right  of  abscissa  indicate  the 
maximum  power  or  amplitude  as  appropriate. 


E- 


CRUSTAL  RESPONSE 


FREQUENCY  (Hi) 

CRUST  IMPULSE 
SP  1 

RESPONSE 

LP 

. ■ . L 

1373 

SEISMOGRAMS 

Figure  15.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  JE-M.  The 
numbers  to  the  right  of  abscissa  indicate  the  maxlmvim  power  or 
amplitude  as  appropriate. 
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Figure  16.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  JR-AZ.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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17.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  KC-MO.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  18.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  KN-UT.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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hikin'  20.  (.riistal  t ri'qiic'ncy  rospoiisc  function,  sliown  as  linear  power  spectrum 
of  the  receiver  crustal  Impulse  response,  receiver  crust  Impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  The 

numbers  to  tlie  right  ol  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  21a.  Crustal  frequency  response  function,  shown  as  linear  power  spec- 
trum of  the  receiver  crustal  Impulse  response,  receiver  crust 
Impulse  responses,  synthetic  seismograms  for  the  LRSM  station 
MN-NV.  The  numbers  to  the  right  of  abscissa  indicate  the  maximum 
power  or  amplitude  as  appropriate. 
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Figure  21b.  Crustal  frequency  response  function,  shown  as  linear  power  spec- 
trum of  the  receiver  crustal  Impulse  response,  receiver  crust 
Impulse  responses,  synthetic  seismograms  for  the  LRSM  station 
MO-ID.  The  numbers  to  the  right  of  abscissa  Indicate  the  maximum 
power  or  amplitude  as  appropriate. 
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Figure  22.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  MV-CL.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  23.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  NP-NT.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  24. 


Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  PC-BC.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  25.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSH  station  SJ-TX.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  27. 
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Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  SW-MA,  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  28,  Crustal  frequency  response  function,  shovm  as  linear  power  spectrum 

of  the  receiver  crustal  impulse  response,  receiver  crust  impulse  » 

responses,  synthetic  seismograms  for  the  LRSM  station  TF—i.L.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  29. 


Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  WH-YK.  The 
numbers  to  the  right  of  abscissa  Indicate  the  maximum  power  or 
amplitude  as  appropriate. 
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Figure  30.  Crustal  frequency  response  function,  shown  as  linear  power  spectrum 
of  the  receiver  crustal  impulse  response,  receiver  crust  impulse 
responses,  synthetic  seismograms  for  the  LRSM  station  WI-NV.  The 
numbers  to  the  right  of  abscissa  indicate  the  maximum  power  or 
amplitude  as  appropriate. 


The  results  of  calculations  show  that  the  total  variation  in  magnitude 
due  to  the  estimated  plane-parallel  crustal  structure  is  much  less  than  the 
observed  variation  in  magnitude  residuals.  The  largest  values  of  crustal 
amplification  are  associated  with  thick  sediments  at  the  surface,  at  Gulf 
Coast  stations  LVLA  and  SJTX  for  example.  Since  tliese  models  are  some 
extreme  it  is  unlikely  that  credible  refinements  in  the  crustal  structures 
would  increase  the  range  of  crustal  amplification.  Major  modifications  of 
structures  can,  in  most  cases,  be  ruled  out  even  on  the  basis  of  the  sketchy 
geological  information  available. 

For  instance,  altliough  addition  of  thin  (d  < 100m),  low-velocity  weath- 
ered layer  or  a thin  layer  of  alluvium  could  be  allowed  at  some  sites  by  the 
available  geological  information,  it  would  be  unlikely  to  introduce  any 
significant  change  in  the  results.  Similarly,  it  would  be  unreasonable  to 
assume  low  velocities  for  sedimentary  rocks  at  some  BUS  stations,  if  it  is 
stated  in  the  site  description  that  these  rocks  are  highly  compact  or  meta- 
morphlzed.  If  a station  located  on  granite  bedrock  the  obvious  choice  is 
to  use  a simple  two  layer  crust  with  a granitic  top  layer.  This  would  be 
incorrect  only  if  the  dimensions  of  the  granite  body  are  extremely  small. 

Tlie  thicknesses  conform  everywhere  with  the  regional  crustal  thicknesses  given 
by  Pakiser  and  Steinhart  (1964),  and  with  the  nearby  crustal  refraction  studies 
and  the  velocities  given  by  Herrin  and  Taggart  (1962). 

Tables  TTa  and  Ilb  summarize  the  results  of  our  calculations.  Table  Ila 

Indicates  the  areal  subdivision  used  (1  for  EUS,  2 for  WS,  3 for  others), 

elves  values  of  density  and  P wave  velocity  at  the  surface  (a  , p ) and  below 

^ s s 

the  Moho  (a  , p ) and  the  raw  vertical  surface  amplitudes  obtained  by  the  Has- 
n n 

kell  matrix  method  for  the  5 kt  and  50  kt  pulses  as  viewed  through  the  LRSM 
o'liort  and  long  period  instruments  (SP5,  SP50,  LP50) . These  include  various 
common  factors  in  the  computations  due  to  pulse  shape,  Fourier  transforms  etc. 
which  are  of  no  Interest  since  only  the  relative  pulse  amplitude  ratios  between 

Pakiser,  L.  C.,  and  J.  S.  Steinhart  (1964).  Explosion  seismology  in  the 
Western  Hemisphere,  Research  in  Geophysics,  MIT  Press. 

Herrin,  E.  T. , and  J.  Taggart  (1962).  Regional  variations  in  Pj^  velocity  and 
their  effect  on  the  location  of  epicenter.  Bull . Seism.  Soc.  Am. , 52 , 
1037-1046. 
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various  stations  matter.  The  magnitude  residuals  of  Booth  et  al.  (l‘)6A)  tor 

sliort  and  long  period  P waves  (SPC  and  LPC)  are  also  given.  Table  1 U>  shows 

reduced  quantities  to  be  analyzed.  We  correct  for  ttie  sub-Moho  medium  by 

dividing  the  calculated  amplitudes  by  a ^ . We  denote  the  quantities  thus 

n n 

corrected  by  tlie  unbarred  designations  SP5,  SP50,  LP50.  Tt»e  base  ten  log- 

arltlims  of  these  are  named  crustal  amplification  terms  in  tills  report.  We 

; — ^ 

name  the  quantity  logj^j^(»'p  ir^)  the  surface  Impedance  term.  To  examine  the 
question  of  whether  a difference  In  observed  magnitude  is  proportional  to  a 
difference  in  the  logarithm  of  theoretical  amplitudes,  the  absolute  level  of 
the  Individuals  comprising  the  pair  of  theoretical  amplitudes  Is  Irrelevant. 

We  are,  therefore,  free  to  subtract  a constant  from  the  logarithms  of  the 
amplitudes  of,  e.g.,  SPS  In  Table  11a,  A different  constant  has  been  sub- 
tracted from  the  logs  of  each  of  the  columns  SP5,  SP50,  and  I,P50,  after  they 
have  been  corrected  for  In  order  to  bring  the  numbers  into  the  range 

of  unity.  The  last  twi>  columns  of  Table  11a  contain  the  magnitude  correction 
terms  given  by  Booth,  Marshall  and  Young  (1974)  for  short  and  long  period 
P-waves.  ISlnce  in  the  following  regression  analysis  we  want  to  test  the 
<lependence  of  the  m<ignitude  correction  on  crustal  structure  together  with  t*; 
in  accordance  with  Her  et  al.  (1977),  we  choose  t*  • 0.2  for  EUS  stations 
anil  t*  “ .AS  for  WHS  stations.  Tliese  averages  are  broadly  valid;  no  reli- 
able t*  figures  are  available  for  individual  stations  at  this  writing.  In- 
dividual station  t*  valiu's  could  he  obtained  by  averaging  many  Individual 
spectral  observations  at  each  station,  and  It  seems  certain  that  they  would 
reduce  the  final  error  variance.) 

In  preparing  for  the  regression  analysis  let  us  first  examine  visually 
some  plots  ot  our  computed  crustal  amplification  terms  vs.  other  relevant 
quantities.  We  show  in  Figures  11  and  32  the  short  period  amplification 
terms  from  Table  lib  (log-amplitude)  plotted  against  the  quantity  logjj^(/p^a^^ , 
where  the  subscripts  s pertain  to  estimated  vales  of  p and  n at  tlie  free 
surface.  The  figures  show  that  the  surface  properties  are  very  important  In 

Der,  2.  A.,  M.  S.  Dawkins,  T.  W.  McKlfresh,  d.  H.  Goncz,  C.  K.  Gray,  M.  S. 

Gllllsple  (1977).  Teleselsmlc  P Wave  Amplitudes  and  Spectra  at  NTS  and 
the  Shoal  Site  as  Gompared  to  Tltose  Observed  In  Eastern  Nortli  America, 
Preliminary  Report,  SDAG-TR-77-9,  Toledyne  Geotech,  Alexandria,  Va. 
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Figure  31.  Short  period  crustal  arapl if icat ion  terra  for  a pulse  from  a 50  kt 
nuclear  explosion  plotted  against  the  surface  acoustic  impedance 
term.  Least  squares  regression  line  is  drawn  through  the  points. 
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Figure  32.  Short  period  crustal  amplification  term  for  a pulse  from  a 5 kt 
nuclear  explosion  plotted  against  the  surface  acoustic  impedance 
term.  Least  squares  regression  line  is  drawn  through  the  points. 
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determining  crustal  amplification  in  our  models.  The  scatter  reflects  the 

deviation  between  amplification  due  solely  to  surface  acoustic  impedances 

and  the  exact  calculation  of  response.  The  scatter  relative  to  a simple 

linear  dependence  is  small.  This  suggests  that  in  most  cases  knowledge  of 

sufficient  to  approximate  the  theoretical  crustal  response.  The 

points  belonging  to  the  EUS  and  WUS  or  third  category  populations  are  mixed 

together  in  these  plots,  showing  that  all  of  these  groups  cover  a similar 

range  in  crustal  amplification  terms,  or  log  (/p  a^.  Since  energy  E •• 

XVJ  s s 

A then  in  a slowly  varying  medium  where  energy  is  concerved  logA  ~ 

-1.0  log/pa  \ Our  calculated  regression  slopes  of  0.327  and  0.393  reflect 
the  fact  that  the  Influence  of  deeper  layers  is  being  felt  on  the  low-fre- 
quency portions  of  the  signal,  and  that  for  low  much  of  t!\e  energy  is 

reflected  before  it  reaches  the  surface. 

As  the  next  step  in  our  analysis  we  tried  to  examine  the  relationship 

SP  LP 

between  the  magnitude  residuals  m^^  and  m^^  of  Booth,  Marshall  and  Young 

(1974)  on  one  hand  and  our  crustal  amplification  terms  log^Pi(SP5),  log^^^ 

(SP50),  log^_(LP50),  log, -(4^  and  t*  on  the  other  hand  for  the  total 
XU  XU  s s 

data  set. 

Figures  33  and  34  show  the  computed  crustal  amplification  terms 
(SP5)  and  logj^^CSPSO)  (called  crustal  amplification  terms)  plotted  against 
the  magnitude  residuals  SPG  of  Booth  et  al.  (1974).  Triangles  denote  stations 
in  western  North  America  (WUS),  squares  denote  stations  in  the  EUS,  circles 
denote  stations  which  could  not  be  classified  into  either  category  (ADIS, 

HWIS  and  WHYK) . Both  figures  show  that  the  points  associated  with  the  EUS 
stations  form  a population  which  is  almost  completely  separated  from  those 

of  the  WUS  stations.  The  magnitude  residuals  are  positive  in  the  EUS  (high 
observed  magnitudes)  and  negative  in  the  WUS.  The  points  for  the  three 
stations  not  belonging  to  either  category  seem  to  group  with  the  EUS  stations. 
Both  the  EUS  and  WUS  populations  exhibit  visible  trends  with  the  crustal 
amplification  terms  although  there  is  a lot  of  scatter. 


The  surface  impedance  term  logj^Q/p^a|  was  also  plotted  against  the 
magnitude  residuals  in  Figure  35.  Since  this  quantity  was  shown  to  be 
highly  correlated  with  the  crustal  amplification  terms  the  distribution  of 
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Clear  explosion  plotted  against  short  period  magnitude  resid 
Is  of  Booth  et  al.  (1974).  Regression  lines  for  the  EUS  and 
S populations  are  drawn  separately.  Regional  magnitude  bias 
shown  on  the  top  of  the  figure. 


lU 


the  points  is  very  similar  to  those  seen  in  the  two  previous  figures. 


In  Figures  33-35  we  also  plotted  the  separate  EUS  and  WUS  lines  derived 
by  regression  between  the  short  period  magnitude  residuals  and  crustal  am- 
plification terms  or  surface  impedance  terms,  respectively.  These  lines  are 
based  on  a regression  analysis  described  later  in  this  report.  The  shift 
between  the  two  lines  represents  an  average  EUS-WUS  magnitude  bias  for  the 
station  networks  used  by  Booth  et  al.  (1974). 

An  anomalous  point  is  MOID,  which  seems  to  group  with  the  EUS  stations 
due  to  its  +0.4  residual  as  given  by  Booth  et  al.  (1974).  This  station  is 
located  in  a missile  silo  on  unconsolidated  silty-clayey  deposits.  In  com- 
puting the  crustal  amplification,  we  assumed  a compressional  velocity  of 
1.5  km/sec  for  the  surface  material.  Although  this  figure  is  quite  low,  it 
could  be  even  lower.  Attempts  to  move  this  point  into  the  VfUS  population  by 
lowering  the  velocity  failed.  Probably  the  difference  reflects  the  effects 
of  the  dense  missile  silo  on  the  silt  "spring".  However,  without  more 
detailed  spectral  analysis  of  this  station  we  cannot  reject  the  hypotheses 
that  beneath  MOID  is  a Q structure  characteristic  of  the  EUS  or  that  a 
powerful  focusing  effect  due  to  non-plane-parallel  layers  is  operative. 


To  derive  the  foregoing  regression  lines  we  set  up  the  following  three 

linear  hypotheses  to  evaluate  the  relative  dependence  of  the  short  period 

SP 

magnitude  resiudals  of  Booth  et  al.  (1974),  on  the  crustal  amplification 

terms  log^Q(SP5),  logj^Q(SP50) , and  quantity  logj^^/p^a^^;  and  on  t*  which  is 
taken  as  .2  for  EUS  and  .45  for  WUS. 


1)  mj^  = a + b logj^Q(SP5)  + ct* 

2)  m^  = a + b logj^Q(SP50)  + ct* 

3)  mp^  = a + b log, a ^ + ct* 

D lU  s s 

The  three  stations  not  belonging  to  either  the  EUS  or  WUS  are  disre- 
garded in  this  analysis.  The  coefficient  of  t*  is  expected  to  be  about 
-iTf/Ln(10)  = -1.36  for  1 Hz  waves.  The  coefficient  is  expected  to  be  1.0 
for  SP5  and  SP50  if  the  calculations  of  crustal  differential  response  are 
performing  on  the  average  as  expected  and  if  the  assumed  source  spectrum  is 
representative  of  the  earthquakes  used  to  determine  njj^  . For  logj^Q/p^o^^  we 
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would  expect  a value  for  b of  roughly  0.3  to  0.4  considering  the  results  of  | 

Figures  31  and  32  where,  for  example,  log(SP50)  = -.327  log/p  a The  ♦ [ 

hypotheses  also  assign  a common  slope  of  b for  both  BUS  and  WUS  populations  in  [ 

each  case.  The  quantity  "a"  is  designed  to  absorb  all  multiplicative  factors,  | 

which  cancel  when  magnitude  differences  are  considered.  ) 

Table  III  shows  the  results  of  statistical  analyses.  The  quantities 
shown  are;  the  linear  regression  coefficients  a,  b,  and  c for  the  three 
hypotheses  with  their  95%  confidence  limits;  the  percentage  of  total  var- 
iance explained  by  each  individual  variable  crustal  amplification  and  t*, 
and  by  their  linear  combination;  and  the  residual  (unexplained)  error  terms. 

The  table  shows  that  the  hypotheses  using  both  variables  can  account 

SP 

for  70-75%  of  the  total  variance  in  magnitude  residuals  Am^^  . The  coefficients 

"b"  are  different  from  zero  at  the  95%  confidence  interval  showing  that 

SP 

crustal  amplification  is  important  in  determining  Am^^  , although  by  itself 
it  accounts  for  only  about  one-tenth  of  the  total  variance.  The  t*  effect 
coefficient  "c"  is  also  significantly  different  from  zero,  accounting  for  ' 

more  than  one-half  of  the  total  variance.  The  coefficient  "c"  ranges  between 
-1.25  and  -1.35  which  is  very  close  to  the  value  -1.36  quoted  above.  | 

The  product  c6t*  is  theoretically  equal  to  [tt  f/Ln(10)  ]6t*.  With  6t*  = 

.25  we  found  for  SP50  c = 1.35  which  implies  f = 1.  Had  we  chosen  6t*  = .17 

(Der  et  al.,  1977)  c would  have  become  1.36  x .25/. 17  = 2.00  which  would  imply 

f = 1.0  X .25/. 17  = 1.47.  Through  March  1978  the  mean  values  of  1/T  at  RKON 

and  0B2NV  were  1.49  Hz  and  1.12  Hz  respectively.  The  mean  frequency  of  1.3 

Hz  is  in  fair  agreement  with  1.47  considering  the  approximate  nature  of  the 

frequency  estimation,  our  lack  of  knowledge  of  the  earthquake  source  spectra 

(which  influence  the  amplitude  and  period  differences  to  a small  degree),  and 

the  approximate  nature  of  the  mathematical  argument  which  is  implicitly  being 

carried  out  in  the  spectral  domain  instead  of  the  more  relevant  time  domain. 

SP 

The  mean  EUS-WUS  Amj^  differential  in  observed  magnitudes  implied  by  our 
values  of  c ranges  between  .34  to  .31  magnitude  units.  The  value  of  the  I 

coefficient  b for  Hypothesis  3 is  smaller  since  the  surface  impedance  term 

varies  over  a wider  range  than  the  crustal  amplification  terms.  It  is  also  ^ 

significant  that  b is  close  to  unity  for  the  50  kt  pulse  which  indicates  that 
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Hypothesis  1 

Independent 
variables  used 


logj^Q(SP5),  t* 

log^Q(SP5) 

t* 


TABLE  III 


Residual  rms 
error  (nij^  units) 

.12 

.20 

.15 


Percentage  of  total 
variance  explained 

70 

13 

53 


Regression  coefficients 
a = -.64 

b = .72  ± .37  (95%  conf.) 
c = -1.25  ± .34(95%  conf.) 


Hypothesis  2 

Independent 
variables  used 

log^Q(SP50),  t* 

logj^Q(SP50) 

t* 


Residual  rms 
error  (m^  units) 

.11 

.20 

.15 


Regression  coefficients 
a = -.83 

b = 1.03  ± .43  (95%  conf.) 
c = -1.35  ± .32  (95%  conf.) 


Hypothesis  3 

Independent 
variables  used 


logio(/^|  ).  t* 
log^o(/pgaj  ) 
t* 

Regression  coefficients 


Residual  rms 
error  (m^^  units) 

.11 

.20 

.15 


Implied  EUS-WUS  regional 
magnitude  bias 
Am^=  . 31 


Percentage  of  total 
variance  explained 

.74 

.10 

.53 


Implied  EUS-WUS  regional 
magnitude  bias 
Am^  = .34 


Percentage  of  total 
variance  explained 

.75 

.12 

.53 


a = .80 

b = .35  ± .14  (95%  conf.) 
c * -1.33  ± .31  (95%  conf.) 
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Implled  EUS-WUS  regional 
magnitude  bias 
Amj^  “ .33 


crustal  amplification  for  this  pulse  is  a proper  estimator  for  magnitude 
residuals.  If  b were  forced  to  a value  of  1.0  in  all  three  hypotheses,  then 
the  highest  proportion  of  variance  would  be  explained  by  Hypothesis  2.  The 
small  value  of  b for  the  5 kt  pulse  shows  that  crustal  amplification  of  this 
pulse  overestimates  the  crustal  effect.  This  is  due,  undoubtedly,  to  the 
high  frequency  content  of  this  pulse. 

SP 

Although  a great  part  (up  to  75%)  of  the  total  variance  in  Am^  can  be 
explained  by  the  combination  of  crustal  amplification  and  t*  in  all  three 
hypotheses  investigated,  the  remaining  variance  is  not  small.  The  remainder 
can  be  caused  by  factors  such  as  the  real  variation  in  t*  in  both  populations, 
imperfect  crustal  modelling,  crustal  focusing,  etc.  It  also  reflects  the 
inherent  variability  of  time-domain  amplitude  measurements  well  demonstrated 
by  variation  of  amplitudes  across  LASA  and  NORSAR,  see  e.g.,  Klappenberger 
(1967),  Blandford  (1974),  Chang  and  von  Seggern  (1977), 

In  Figure  36  we  see  that  the  observed  variation  of  long-period  residuals 
is  much  greater  than  can  be  explained  by  crustal  structure,  and  there  is  no 
EUS-VTUS  separation.  Evidently  focusing  and  other  effects  are  dominant. 


Klappenberger,  F.  (1967).  Distribution  of  short-period  P phases  over  LASA, 
Seismic  Data  Laboratory  Report  287,  Teledyne  Geotech,  AD  815-580. 

Blandford,  R.  (1974).  Short  period  signal  to  noise  ratio  at  NORSAR,  SDAC-TR- 
74-13,  Tfeledyne  Geotech,  ADA  015650. 

Chang,  A.  and  D.von  Seggern  (1977).  A study  of  amplitude  anomaly  and  m^  bias 
at  LASA  subarrays,  SDAC-TR-77-11,  Teledyne  Geotech,  Alexandria,  Va. 
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